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We analyse the time of equilibration of subsystems of a larger system under a random total 
Hamiltonian, in which the basis of the Hamiltonian is drawn from the Haar measure. We obtain 
that the time of equilibration is of the order of the inverse of the arithmetic average of the Bohr 
frequencies. To compute the average over random basis, we compute the inverse of a matrix of 
overlaps of operators which permute four systems. We first obtain results on such matrix for a 
representation of an arbitrary finite group and then apply it to the particular representation of the 
permutation group under consideration. 



I. INTRODUCTION 

The phenomenon of convergence to equilibrium despite 
an underlying deterministic dynamics was usually justi- 
fied by referring to subjective lack of knowledge, i.e. by 
putting probabilities by hand. However, already in 1929, 
von Neumann (see [l[ for an English translation and com- 
mentary) put forward an argument for relaxation with- 
out referring to an ensemble: for a typical initial pure 
quantum state, averages of macroscopic observables will 
be for most time around their equilibrium value. In this 
approach, thermalization is implied by statistical prop- 
erties of quantum states themselves, namely it is due to 
the fundamental lack of knowledge represented by quan- 
tum probability. This "individualist" approach to equi- 
librium (as phrased in 111) has been recently intensively 
developed, see e.g. (JtSI- More broadly, new theorerical 
and experimental developments in the question of sub- 
system equilibration in close quantum systems have also 
been achieved However the time of equilibration, 

a very important aspect of equilibration and thermaliza- 
tion, was not considered so far. A natural time scale that 
appears from the analysis of [3| is the inverse of the small- 
est energy gap of the Hamiltonian. However the latter 
is typically exponentially small in the size of the system, 
and thus cannot offer an explanation for the fast nature 
of thermalization. 

In this paper we consider the issue of time of equili- 
bration. As in [i| we consider system S and the bath 
B, and we are interested in equilibration of the system, 
given that the bath is sufficiently large. We evaluate the 
distance of the state p(t)sB of a system and a bath, evolv- 
ing according to a random Hamiltonian, from the state lo 
which is obtained by removing the blocks of p(0) which 
are off-diagonal with respect to Hamiltonian spectral de- 
composition. 

Our main result amounts to showing that if we choose 
the basis of the Hamiltonian randomly according to Haar 



measure, then the time of equilibration depends on the 
(arithmetic) average distance between the eigen-energies 
of the Hamiltonian rather than on the worst case gap. 

Computing the average over random choice of the 
eigenbasis of the Hamiltonian, is reduced to evaluating 



average of the sort Tr 



U m X{U®yY 



over Haar dis- 



tributed unitary transformations U, where X, Y are op- 
erators. This leads us to a general problem of inverting 
a matrix M g h = x(<? -1 ^)' where g, h are elements of a 
finite group G and \ i s character of some given repre- 
sentation of the group. It turns out that such a matrix 
enjoys nice properties, which allow to obtain the inverse 
in the case of interest (i.e. for G = S4). We also present 
some other properties of the above matrix. 

The paper is organized as follows. In sec. [II] we com- 
pute average of the distance from equilibrium over ran- 
dom basis of the Hamiltonian. In sec. Mil we discuss the 
dependence of time of equilibration on eigenvalues of a 
Hamiltonian with a random basis. In sec. IIVI we derive 
the inverse of the matrix of overlaps of swap operators. 
To this end we first provide properties of a general matrix 
of this type constructed from any given representation of 
a group. 



II. AVERAGING OVER RANDOM CHOICE OF 
EIGENBASIS 



Let us introduce some notation. We have two sys- 
tems S and B (the latter one playing the role of the heat 
bath). We consider an arbitrary initial state psb(O) = 
|V ; )s-b(V'I- Since we shall consider random Hamiltonians, 
whose basis is chosen according to Haar distribution, we 
can equally well take a standard product initial state: 
IV'sb) = |0)s|0)b. We now consider the evolved state 
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Pss{t) given by 

PSB(t) = e- lHt p SB (Oy m 
We also define a state u>sb as 

i 

where Pi are eigen-projectors of the Hamiltonian 



We set 



H = UH U\ P t = UP°U\ 



(1) 
(2) 

(3) 
(4) 



where Hq = diag(E\, . . . , Ed). We assume that the prob- 
ability measure over Hamiltonians splits into two parts 



d/i( J ff)=d/x 1 (C/)d/i 2 ( J ffo) 



(- r » 



where p\ is Haar measure, while p2 is distribution over 
energies (such a separation holds e.g. for gaussian unitary 
ensemble). Therefore we have (. . = (((. . .)u))h Let 



us also denote W 



AHt 



Wo 



oil H a 



W ® W T , 



and P = Ei-Ff ® -Pf- B Y Vx i: x 2 we will denote the 
operator which swaps the systems X\ and -X" 2 . 

We consider the distance between the state ps(t) = 
TiBPBs(t) and the state u>s = Ttb{u>sb), induced by 
the Hilbert-Schmidt norm \\A\\ 2 = y/Tr(A^A), and we 
want to average it over random Hamiltonians. We will 
need the following proposition 

Proposition 1 The following relation holds 



Tr 



YU m XU^ 



(6) 



where the label i denotes the system SiBi, which are iso- 
morphic to SB, X is given by 



X = (W-P)i 3 ® (W f -P) 24 
with a = Psb(0), o-nd Y is given by 

Y 123 4 = Vi2:34(CTl ® 0- 2 ® F34). 



(7) 

(8) 



To prove this proposition, we shall need the following 
easy to check fact 

Fact 2 For two systems 1 and 2, and for arbitrary oper- 
ators Ax, B2, C12 and D12 we have 

TrlCiaAi ® B 2 D 12 A\ ® Pj] = 

Tr[V 12 :34(Ci 2 ® L> 34 )(Ai ® P 2 ® 4 ® b\)] (9) 

i/ere 3 and 4 are auxiliary systems isomorphic to 1 and 
2. 



Proof of Proposition [TJ We rewrite Tr{ps{t) — W5) 2 
as follows (we will not write dependence on time explic- 
itly) 

Tr(p s - u> 3 ) 2 = Trp| - 2Trp s cj s + Tru 2 s = 

Tr [(PSi ® Ps 2 - P5i ® w S2 - W Sl ® p5 2 + UJ Sl ® UJ Sl )V Sl :S 2 ] 
= Tr( Pl ® p 2 F) - Tr( Pl ® u; 2 F) - 

-Tr(wi ® p 2 F) + Tr(wi ® w 2 F) (10) 

where Fi 2 = Vs t : s 2 ® s 2 , an d again, the label i denotes 
system S l B i , so that e.g. pi = ps 1 B 1 - 

Consider now the first term. Writing psb = 
e -iHt asBe iHt we obtain 



,-ifft iHt 



Tr(pi ® p 2 ¥) = Tr(e ""ae 
= Tr(o-i ® cr 2 Vyi ® W 2 ¥W} ® W 2 + ) 



e - iK Ve ii? *F) 



fill 



We can now use Fact [2 putting Ci 2 = crj ® cr 2 , Z?i 2 = F, 
A = B = W . As a result we get 



Tr(pi ® p 2 F) = Tr[Fi234(^i ® W 2 ® W x f ® W 2 + )] 



= Tr 



ll234^® 4 (W / ® ^0 ® ® W,})!/* 



®4 



(12) 



In similar way we get 
Tr(pi ® w 2 F) = 



r i2 34C/® 4 (Ty ® Pi ® W f ® P;)^ 



Tr(wi ® p 2 F) = 



n 23 4f/® 4 (P ® Wo ® Pi ® Wtf )tf +S 



Tr(wi ® w 2 F) = 



Y 127ii U®\P l ® P; ® Pi ® P,-)C/ tS 



(13) 



If we now insert (fT3")) and (|T2"T) into (|10[) we obtain ([5]). | 
We are now prepared to state the main result of this 
section 



Theorem 3 We have up to a leading order 



M 2 1 



d 2 d s 



Ps ( t )- Ws || 2 V~^ r - + (^-^) +o{-) 



d 2 d 2 



1 



whe 



(14) 



£ = Wo = E l die^\ V = Triy 2 = ^ d. v 
7 = Y,i d h d = d sdB 



T2E,l 



(15) 



and (•)[/ denotes average over U according to Haar mea- 
sure, and di are the dimension of the degeneracy. 

Remark. It is easy to show that in the case of no- 
degeneracy, average from Theorem [3] reduces to 



\p s (t) - u) S \\l)u 



1pTs + ^ + 0{ T b 1 (16) 



3 



where in this case £ = TrW = Tre tHot and 77 = TrW, 2 = 
Tre i2Hot . 

Proof of the Theorem \S[ We want to compute 
the average over U of the formula ©• Let us denote 

t 4 (-) = / dUU m {-)U m] . We thus want to compute 
Tr[YY 4 (A)] with X,Y given by and © respectively. 
To this end we apply Fact [23] from the Appendix. We 
write our operators X and Y in terms of product of op- 
erators as follows: 

X = C\-C 2 -C a + C 4 , (17) 

where 

C 2 = w o ® Wo ® p ° ® 

i 

C 4 = ^2P?®P°® P° ® P? (18) 

where we put systems in the following order (3, 4, 1, 2). 
We also introduce Y as follows 

y' = °i ® ^ ® 4' ® 4* (19) 

Now. we define the vectors c^ fe ^ by = TrC^V^-i, 
fc = 1, . . . ,4 and a by a„ = Tr(yV„.-i). Note that, we 
present here, result for a a which is related to the vector 
a, in a similar way, like Y to Y , namely a = Vi2 : 34a . In 
order to compute these coefficients we decompose a given 
permutation into cycles, so that for product operators, 
as those in the above equations, the formula breaks into 
product of separate terms associated with the cycles. For 
a single cycle we use Fact [22] We then obtain 

ci = m\ \ed, \ed, \e, ici 2 ,^ 2 , i^, d\ d, d, i^i 2 , 

iei 2 ,d^f,iciM^i 2 ,d,d,i^i 2 ,ieiMei 2 d,d,d 2 ), 
c 2 = ( 7 i£i 2 ,d,pr,7,7,^,dici 2 ^ 2 ,iei 2 ,d,d,iei 2 , 

|£| 2 ,d, P r,7,7,d,d,l^| 2 ,7,P*e,rf,7), 

c 3 = (7iei 2 ,rf,K*,7,7,p*e,rfici 2 ,d 2 ,iei 2 ,d,d,iei 2 ,i£i 2 ,d, 

PC 7> 7, d , d, |£| 2 , 7,p*£, d, 7), 
c 4 = (7 2 , 7^, t, 7, 7, t, 7d, d 2 , 7, d, d, 7, 
7, d, t,7,7,d, d,7,7,i, d,7), 

a = (ddB,dds,dB,ds,ds,dB,ddB,dds,dB,ds,ds,dB, 

d B ,ds,d B , d s , 1, 1, ds,ds,ds,dB, 1, 1). (20) 

Here, $ = Tr(pWo) = d 4 e^', £ = TrW = £<6 = 
'^2 i die lEit , p = di£i where dj are the dimension of 
the degeneracy, and finally "the scalar product": p = 
= Ei^fe^'. What is more, 7 = £\ d? and 

Now, we proceed to compute the matrix M l . We 
refer to the section IIVI where we consider the general 



properties of the matrix M defined more generally for 
representation of any group. In our case, for d > 4 the 
matrix is invertible, and its inverse is given by (|66p . We 
can now use the formula (|23p and obtain 

(\\ Ps (t) - US ^u = - da(2 + d)( _ 3 ! d + 3( p + ( p ) >< 

(id + Ad 2 + 2d 3 - 4dd B + M 3 d B + d i d B - 
|£| 4 (2 + d)(l + d - ds - d s ) - Adds - 2d 3 d s - 
4d 4 d s - d 5 d s - 6(2 + (-2 + 2d + 4d 2 + d 3 )d B - 
(2 + 4d + d 2 )ds) + 4 7 - 4d 7 - 6d 2 7 - 2d 3 7 - 
4d B 7 + 2dd B 7 + 5d 2 d B 7 + d 3 d B 7 - 4d S 7 + 
2dd S 7 + 5d 2 d s7 + d 3 d S 7 - 2 7 2 - 3d7 2 - 
d 2 7 2 + 2d B7 2 + dd S 7 2 + 2d s7 2 + dd s7 2 + 
2|£| 2 (1 + d-d B - d s )(d + 2 7 + dry) + 

—2 —2 —2 

4dt — 4d^t — 4dgt + £ rj + d£ 77 — d^£ — 

dsfv + ?V + dfrj - d B e 2 77 - d s £ 2 r] - 4p*£ - 
4dp*£ + 4d B p*£ + 4d s p*£ - 4p£* - 4dp£* + 
4d BP r+4d sP r) (21) 

One then finds, that up to a leading order, and a constant 
factor, this gives the right hand side of (|T4|) . | 

Remark. Using ideas of measure concentration it is 
easy to show that Theorem [3] can be extended to say that 
the vast majority of unitaries U will have a — Ws|| 2 

close to the average and hence the correspondent Hamil- 
tonian will equilibrate quickly. Moreover, we can pass 
to trace norm, by using ||^4||i < \/I}||A||2, for operators 
A acting on C D , which adds factor ds (recall that we 
consider ds ~S> ds). 

At the end of this section, let us present important re- 
mark, but first we will need the following lemma, namely 
Levy's lemma [l8| 

Lemma 4 For a given function f which is Lipschitz's 
continuous with a Lipschitz's constant, it can be shown 
that: 

Pr (\f(U)-v\>S)<Ce~ cds2 , (22) 

where C, c are constants and d is the dimension of a total 
system. 

We expect exponent from equation (f2"2"j) to be small, so 
we put S = d~3 in our case. 

Remark Application of Levy's lemma and passing to 
trace norm from the above remark yields that for our 
average from Theorem [3J we have that with high proba- 
bility, according to Haar measure over U, it holds that 

\\ps(t) - ws||i < 

where c is an absolute constant, and the other notation 
is as in Theorem [3J 
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III. AVERAGE OVER ENERGIES 

Having obtained the expression which depends only 
on eigenvalues let us first observe what happens if we 
average over time. Then using Eq. (|14p and averaging 
over a fixed time interval [0, T] we find 



< 



1 

1 



dt{\\p s (t) - u s \\l)u 

\ e i(E 3 -E k )T _ -H 



(24) 



TdPds ~f 

3,k 



Td 4 

j,k,r,a 



\Ej-E k \ 

,i(Ej —E k + E r — E a )T _ j | 

\Ej - E k + E r - E a \ 



Since 



and likewise 



: i(Ej-E k )T _ x | 
\Ej~E k \ 



< 1 



AE 3 -E k + E T -E S )T _ j| 

\Ej — Ek + E r — E s \ 



< 1 



(25) 



(26) 



we see that if we choose T of order of the average (over 
the indices j, k, r, s) of the inverse gaps (Ej — Ek)' 1 and 



(Ej - E k + E r - Es)- 1 , the LHS of 423} will be small. 
This in turn ensures that the state of the subsystem is 
close to los for most times in the interval [0,T]. 



A. Average over gaussian independent energies 

As a simple example, let us now consider the average 
over Gaussian independent energies, i.e. we consider the 
following distribution 



where 



e {E 1 ,,..E d ) = g(Ei) . . . g(Ed) 



Q{E l ) 



(27) 



(28) 



and a = log d. We obtain 



Proposition 5 For the ensemble of Hamiltonians, ob- 
tained by picking a basis according to Haar measure and 
energies according to distribution (|27p we have 



e -2t 2 (logd) 2 + 0(_L) < (||p 5 (t) - U S \\1)b 



< 



d B 

-4t 2 (io g d) 2 +0 r_Ly 
ds 



(29) 



Remark: It follows that the time of convergence of 
the state p{t) to equilibrium scales roughly as the inverse 
of log d, i.e. as the inverse of the volume of the total sys- 
tem. 



Proof: Due to Theorem [3] in order to prove the propo- 
sition we need to compute the average (from now on, we 
omit the lower index in averages) 



d 2 d s 



d 2 d 2 



(30) 



over the distribution (|27|l . Since for any random variable 
Z we have {Z 2 ) > (Z) 2 , we obtain 



d 2 



< 



< 



\W 
d 2 



\vf_J_ 

d 2 d s 

nil 
d 4 



d 2 



d 2 



< 



(hi 4 > , m 4 



f / 4 



(31) 



since we can drop term i-^s) , because it is a constant in 
our case and 7 2 |£| 2 ~ \/I?F- We first compute (|£| 2 ): 



JE m t„-iEnt 



E m t 



He 



-iE n t\ 



where we used independence of energies. Since (e lEnt ) 

4-1 2 

e 2 we obtain 
(Ifl 2 ) 1 



(32) 



d 2 



d 2 



(d + (d 2 - d)e~ 2t2 ° 2 ~ e- 2tV2 + 0(h (33) 

d 

Setting a = log d an dropping O ) we obtain the lower 
bound. To get the upper bound we need to perform a 
similar, but a bit more tedious, calculation which gives 
the following expression for (|£| 4 ) (the full calculations 
are a little bit too lengthy and rather uninformative, so 
they are not reported here) 



ML 4 ) v 1 
d 4 ^ d 4 

klmn 



e iE k t e -iEit e iE m t e -iE n t\ 



E 



(e lEkt ) (e- lElt ) (e lEmt ) (e"^"*) + O(-) 



= e 



-Wo* 



(34) 



Now, since 0(±) < 0(-^), and r}(t) = £(2t) we obtain 
the upper bound. 



IV. INVERSE OF THE MATRIX M 

In this section of derive properties of the matrix M 
which were needed in the proof of Theorem ([3]) . 



A. Properties of M matrix for general 
representations 

We will first introduce some notation. Denote by G an 
arbitrary finite group, \G\ = n. Let 

D a : G -> Hom.(n a ); a = 1, 2, r; dim?T = d a (35) 
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be all inequivalent, irreducible representations (IR) (not 
necessarily unitary) of G and let 

D a (g) = (DgG/)); i,j = 1, 2, d a (36) 

be their matrix forms where D 1 (g) = 1 is the trivial 
representation. By 



X a {g) = Tv{DUg)) 



(37) 



we denote corresponding irreducible character (ICH). 
We now define our main object - the matrix M D . 

Definition 6 Let D : G — > ffom(K) be any represen- 
tation (not necessarily unitary) of G. Define a matrix 
M e M(n,C) 

M D = (m gh ) = (Tr(D-\g)D(h))) 

= (MDig^h))) = (x^" 1 /*)) (38) 

We apply this definition to irreducible representations 
D a : 

Definition 7 For irreducible representations D a we de- 
fine the corresponding matrices 

M a = KJ = (Tr(D a )-\g)D a (h))) 

= (Tr(D<*(g- 1 h))) = ( X a (g- 1 h)) (39) 

Thus from the definition of M a , it follows that in order 
to calculate the entries of M a we do not need to know 
explicitly irrep D a , but only ICH \ a ■ 

Now we shall express the matrix M D by means of the 
matrices M a . Namely, from the decompositions 



D = ® r a=1 k a D a ; k a eNU{0} 



n 



X~ = k a x" 

a=l 

(40) 

where k a is the multiplicity of irrep D a in D and from 
the character properties we get 

Proposition 8 1. Matrices M a are Hermitian and 

r 

M D = J2 k *M a =*> (M D ) + = M D . (41) 

a=l 

2. The sum of elements in each row and column of the 
matrix M D is equal to nk\. 

Further, using orthogonality relations for ICH 

-Ef^r 1 )^^ (42) 
which one can derive from Schur's lemma, one can prove 



Proposition 9 The matrices M a are proportional to or- 
thogonal projectors: 



M a M> 3 = —5 al3 M c 

dry 



(43) 



whereas the matrices P a = ^-M a form the complete set 
of orthogonal projectors: 

r 

pap/3 = gafipa. P a = 1; (P a ) + = P a (44) 

a=l 

In particular the matrices M a and P a mutually com- 
mute. 

This already gives us eigenvalues of the matrix M D in 
terms of dimensions d a and multiplicities of the irreps, 
which allows to derive the formula for the inverse of Mb , 
whenever it exists (see Theorem [T5"|) . We can however 
also find eigenvectors in terms of matrix elements of ir- 
reps. Namely, consider n vectors in C™ whose entries 
are defined by the matrix elements of irrep D a in the 
following way 



U§ = (D^g- 1 )) G C": 9 eG; a = 1,2, ....,r; 
i , j ' = 1 , 2 , — , d a . 



(45) 



where a, i, j label the vectors and g £ G label the 
entries of the vector [/"■ <E C™ i.e. the vector [/"■ has the 
form 



{U^ 1 = {Dt j {gi%Dt j {g^),... i Df j {g- 1 )) e C" 



and in particular 



(f/ 1 ) T = (i,l ) ...,i)e 



(46) 



(47) 



It turns out that these vectors are eigenvectors of the 
matrices M a : 

Proposition 10 The are linearly independent and 
they are eigenvectors for matrices M a and P a i. e. 

M a U% = S a0 U% ; P a U% = 8 af) 17% (48) 

If the irrep D a are unitary then the vectors U-j are or- 
thogonal with respect to the standard scalar product in 



Proof: We prove the proposition in Appendix. 

As an easy corollary from Proposition [10] we get the 
following theorem concerning the eigenproblem for the 
matrix M D . 



Theorem 11 The vectors are eigenvectors for the 

matrix M D i.e. 



(49) 



6 



and the eigenvalues of M are the following 

\p = kf) — . (50) 
dp 

The spectral decomposition of M D thus reads 

r 
a=l 

where the eigenprojectors P a are defined in Proposition 
® 

Directly from this theorem it follows 

Corollary 12 1. The matrix M D is invertible iff each 
multiplicity k a in the decomposition 

r r 

X D = ^k aX a M D = Y J k a M a (52) 

a—1 a— 1 

is nonzero. 

2. For a given a the vectors U?- i,j = 1,2, ....,d a 
span the eigenspace for the eigenvalue X a , so the multi- 
plicity of A Q is equal to d a . 

3. The eigenvectors U°j does not depend on the rep- 
resentation D : G — > Hom(V), whereas the eigenvalues 
X a depends on the representation D : G — > Hom(V) via 
multiplicities k a . 

4- We have also 

detM D = U r _(k a ^-) dl , 

d a 

r 

TtM d = ^""^ nkgd a = n dim D 

a = l 

Thus in order to calculate the eigenvalues A a of the 
matrix M D we need only the multiplicities k a of irrep 
D a in the representation D (the dimensions d a and rank 
n = \G\ are known). 

From the above spectral decomposition we get 

Corollary 13 // the matrix M D = J2 a =i k a M a is in- 
vertible k a > 1) then 

a— 1 a — 1 a—1 

(53) 

In fact this formula express the entries of the matrix 
(M ) -1 in terms of ICH, namely we have: 

i M °)~ g h = ^Er^Gr 1 ^ (54) 

a — 1 

i.e. all we need to calculate (M 11 )" 1 are ICH and the 
multiplicities k a of irrep D a in the representation D. 



Remark 14 It is known fldl] that one can calculate the 
multiplicities k a of irrep D a in an arbitrary representa- 
tion R of the group G using the following formula 

k a = (x R , x a ) = - E x R (g)x a (g^), (55) 

see 

where (x R iX a ) * s the scalar product in the linear space 
of central functions on the group G. 

Finally, we want to express the inverse of M D as a 
polynomial of M D . To this end, note that from the Her- 
miticity of the matrix M D it follows that the rank of the 
minimal polynomial of M D is equal to r and the coef- 
ficients of this polynomial are determined by r pairwise 
distincts eigenvalues of M D . Thus it is possible to write 
the matrix (M )" 1 as a polynomial of degree r — 1 in 
M D . In fact we have 

Theorem 15 Let 

W(x) = x r + s r _ix r_1 + ... + sxx + sq (56) 

be a minimal polynomial of the matrix M D i.e. 
W(M D ) = 0. Then if s ^ 

{M D )- X = —[(M D y- 1 + Sr _ 1 (AI D ) r - 2 + ...+s 2 M D +s 1 } 
so 

(57) 

This formula express the inverse of the matrix M D as a 
polynomial function of itself. 

In next section we shall apply these results to our repre- 
sentation. 

B. Applications 

In this subsection we will apply the above results to 
a particular representation of the symmetric group S n 
which is described in 

Definition 16 Let H = <8>™ =1 C d , so dimH = d n . We 

define the representation D of the group S n in the space 
H by mean of operators which swap subsystems: 

Vcr e S n D (o-)(e h <g> e, 2 ® ... ® e» n ) 

= ® e ff -i( i2) ® ... ® e CT -i (in) (58) 

where {ei}f =1 is a basis ofC d . 

In other words D(a) = V a , using notation from previous 
sections. 

An important property of any representation is its 
character and in this case it is not very difficult to prove 
that 

Proposition 17 The character of the representation D : 
S n — > Hom(H) has the following form 

WaeS n x D (a)=d l( ~°\ (59) 

where 1(a) is the number of cycles in the cycle decompo- 
sition of a G S n . 
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It follows that in the case of the representation D of 
S„ the matrix M D has the form 



M 



D 



Example 18 For the group S3 the matrix M 
following 



(60) 
is the 



M 



D 



Id 3 d 2 d 2 

d 2 d 3 d 

d 2 d d 3 

d 2 d d 

d d 2 d 2 

\ d d 2 d 2 



d 2 
d 
d 
d 3 
d 2 
d 2 



d 
d 2 
d 2 
d 2 
d 3 
d 



d \ 
d 2 X 

d 2 
d 2 
d 

d 3 ) 



(61) 



From Theorem [TT] and Corollary [T2] of the previous 
subsection it follows that in order to describe the basic 
properties of the matrix M D , in particular its eigenvalues 
and the inverse (M D )~ l , one has to calculate the multi- 
plicities k a of irrep D a in the representation D. Using 
the formula from Remark Q3] and the character tables for 
S 3 and 5*4 [Hi one gets 

Proposition 19 1. The multiplicity coefficients k a for 
S3 are the following 

h = -(d 3 +3d 2 +2d); k 2 = -(d 3 -3d 2 +2d); k 3 = -(d 3 -d) 
6 6 3 

(62) 

2. The multiplicity coefficients k a in the case of S4 
are of the form 



h = 


^d(d + 


l)(d + 


2)(d + 


3); 


k 2 = 


\,d{d- 


l)(d- 


2)(d- 


3); 


k 3 = 


\d\d 2 


-i); 






k± = 


\d { d 2 - 


-i)(d 


-2); 




h = 


\d { d 2 - 


-i)(d 


f2) 





(63) 

From Theorem[TT]we get immediately the values of the 
corresponding eigenvalues and then from the Corollarv[T2l 
and Theorem [TS] we get 

Theorem 20 1. For S3 we have 
1 



where 
ail = d 6 



CI14 
021 
024 
031 
034 
O41 
O44 
O51 
O54 
061 
064 



d 3 - 
d 3 - 

2d 4 
d 3 - 
2d 4 - 
d 3 - 
d 6 - 
2d 4 - 
d 3 - 
2d 4 
d 3 - 



3d 4 + 2d 2 ; 012 
d 5 ; 015 = 2d 4 - 
d 5 ; o 22 = d 6 - 
-2d 2 ; a 25 = d 3 
d 5 ; a 32 = 2d 4 - 
-2d 2 ; a 35 = d 3 

- d 5 ; a 42 = 2d 4 - 
-3d 4 + 2d 2 ; a 45 
-2d 2 ; a 52 = d 3 

- d 5 ; a 55 = d 6 - 
-2d 2 ; a 62 = d 3 

- d 5 ; a 65 = 2d 4 - 



= d 3 
-2d 2 
3d 4 - 
-d 5 
-2d 2 
-d 5 
-2d 2 
= d 3 
-d 5 
3d 4 - 
-d 5 
-2d 2 



-d 5 ; 
ai6 
2d 2 ; 

026 
033 
036 
043 

-d 5 ; 

053 
2d 2 ; 

063 
a 66 



Ol 3 = 
= 2d 4 

023 = 
= d 3 - 
= d 6 - 
= d 3 - 
= 2d 4 

O46 = 
= d 3 - 

O56 = 
= d 3 - 
= d 6 - 



d 3 -d 5 ; 
-2d 2 ; 
2d 4 - 2d 2 ; 
d 5 ; 

3d 4 + 2d 2 ; 
d 5 ; 

-2d 2 ; 
d 3 -d 5 ; 
d 5 ; 

2d 4 - 2d 2 ; 
d 5 ; 

3d 4 + 2d 2 



and s 3 = d 3 (d 2 - l) 2 (d 2 - 4) = 9d 5 - 4d 3 - 6d 7 + d 9 
Theorem 21 2. For S4 we have 



M~ l = —(A/ 4 - siM 3 + s 2 M 2 - S3M 1 + s 4 l) (66) 
S5 



where 








81 


= d 2 (5d 2 


+ 19); 




S2 


= 2d 2 (d 2 


- l)(5d 4 + 23d 2 - 


K20); 


S3 


= 2d 4 (d 2 


-l) 2 (5d 4 + 7d 2 - 


1-12); 


S4 


= d 4 (d 2 - 


l) 3 (d 2 -4)(5d 4 


-9d 2 + 36); 


S5 


= d 6 (d 2 - 


l) 4 (d 2 -4) 2 (d 2 


-9) 



(67) 



and in this case d ^ 1,2, 3. 



M- 1 = 



d 3 (d 2 - l) 2 (d 2 -4) 
where d 7^ 1 , 2 and 



(M 2 -3d(d 2 +l)M+3d 4 (d 2 -l)l) 
(64) 



1 



Ar 1 = — x 

S3 

I an 012 ai3 014 015 aie ^ 

a 2 i a 22 a 2 3 a 2 4 a 2 5 a 2 6 

031 03 2 033 034 035 036 

041 a4 2 043 044 045 046 

051 052 053 054 055 «56 

V 061 «62 063 064 065 066 / 



(65) 



V. CONCLUSIONS 

We have shown that the time of equilibration of a small 
subsystem under the dynamics of a random Hamiltonian 
is fast, being determined by the mean inverse of the 
energy gaps of the Hamiltonian, which in typical cases 
scales as the number of particles in the system. This 
should be contrasted with the time scale that can be ob- 
tained from the results of [|J , which is given by the inverse 
of the smallest energy gap of the Hamiltonian. The main 
message of this work is that in order to understand the 
time of equilibration in quantum systems, one must con- 
sider more than the eigenvalues of the Hamiltonian. In- 
deed, the structure of the eigevectors of the Hamiltonian 
appears to be of crucial importance for equilibration to 
happen quickly. Interestingly, asymptotic equilibration 
can be inferred just from the knowledge of the eigenval- 
ues of the model, this being the main result of [J]. 

In our work we have shown that for almost any choice 
of the eigenvectors (when picked from the Haar mea- 
sure), equilibration will happen quickly. A direct conse- 
quence of our result is that we can replace the Haar mea- 
sure when choosing the basis by any quantum unitary 4- 
design, since we only used averages over four moments of 
the distribution in our arguments. As random quantum 
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circuits of order n 4 gates form a unitary 4-design (20j . 
this means in particular that most Hamiltonians whose 
eigenbasis are determined by a sufficiently large quantum 
circuit (with more than 0(n 4 ) gates) are such that small 
subsystems equilibrate fast. A drawback of the result is 
that typically a Hamiltonian chosen in this way will be 
very different from realistic Hamiltonians, which should 
be formed by a sum of few-body terms. It is an inter- 
esting open problem if one can say something about the 
generic case of this more realistic type of model. 

Note added: After the completion of this paper we be- 
came aware that similar results have been reported in 
[H| and [13. 
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VI. APPENDIX 

We prove here a few auxiliary lemmas and facts. 
Fact 22 For ir G S n being a cycle, we have 

Tt{V n A x ® . . . g> A n ) = Tr{A n(1) . . . A <n) ) (68) 
Proof. By direct inspection. 

Fact 23 Consider the twirling operation r„ given by 
t„(-) = / dUU^ n (-)U^ . Then for any operators A and 
B acting on (C d )® n we have 



Lemma 24 Let {tpi} be arbitrary set of vectors from the 
Hilbert space Ji. Let M be the matrix of the elements 
from the set: Mij = (tpi\ipj), and denote by M^ 1 the 
pseudoinverse of M , i.e. the unique matrix satisfying 
M M = MM = Q, where Q is an orthogonal projec- 
tion onto a support of the matrix M (Q is the orthogonal 
projection onto the range of M). Then the orthogonal 
projector P onto the subspace spanned by {ipi} can be 
written as: 



(71) 



where Xij are elements of matrix X and by X we mean 
X = M , so the pseudoinverse of matrix M. 

Proof. 

We must show that the operator P is indeed an orthog- 
onal projection i.e. that P = P 2 . Let start our proof by 
writing the following expression for the P 2 : 



J2x tJ x kl \^)(^ k )(in\ 



ijkl 



ijkl 



= J2x ij \^ l \J2M jk X kl , (72) 

ijl k 

where we use definition of M from Lemma We can 
now express our equation in terms of Q and use this to 
obtain desired result 

p 2 =Y j x 1] \^)(MQ 3 i = E (»< I WiX^I 

ijl il \ j I 

= Y,(XQ) u \i>i)(ipt\ =Y,Xu\i>i)(M = P, (73) 

il il 

since according to Lemma [5U MQ = QM = M and 
XQ = QM = X. | 



TrL4T„0B)] = (alM-^b) 



(69) 



where a = (a 7r ) neSn , b = (M^eS™; with a ^ = 
TtAV^-i jb^ = TtBV w -i. The matrix M is given by 
M W , CT = (V K \V tr )=Tr(y n -iV a ). 

Proof. It is easy to check that the twirling operation 
is an orthogonal projector in the Hilbert-Schmidt space 
of operators, with the scalar product (A\B) = Ti(A'B). 
It projects onto the space spanned by the permutation 
operators V w . Then from Lemma [2"4l we have that 



Tt[A*t(B)] - ^2,{A\V 7! ){M~ l ) Vi(T {V„\B) 



(70) 



However (A\V n ) = Tr(AW v ) = a* and similarly 
(V a \B) = b fj , where * stands for complex conjugate. This 
ends the proof. 



A. Proof of proposition 1101 

In order to prove this Proposition we will need: 

Lemma 25 Let \ : G — > C be any character of the group 
G (or even any central function on G) and D a be an irrep 
ofG. Then 



gee da 



where (•, •) is a scalar product in the space C c 

Now we can prove the Proposition 1101 
Proof. 

(M a U^ g = J2x a (9^h)D^(h-') 



(74) 



(75) 
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we set 

u' 1 =g- 1 h (76) 

then 

(M°Uf>) B = E^K^X-C^ 1 ) (77) 

u 

now we use the above Lemma and the fact that ICH of 
G are orthonormal i.e. (x a iX^) = an d we get 

= (78) 

B. Miscellaneous facts about matrix M D 

It turns out that the matrix M D may be written as a 
linear combination of adjency matrices of the so called 
Commutative Association Scheme (see (23|) determined 
by the class structure of the group G. 

Definition 26 Let C\ = {e}, Ci t ....,C r be the conjugacy 
classes of the group G. We define the i th relation Ri on 
G x G in the following way 

(g,h)eRi ^ g~ 1 heC l . 



Then the pair (G, {i?i}[ =1 ) is a Commutative Assotiation 
Scheme and by A4 we denote the coresponding adjency 
matrices which are matrices of degree \G\ = n whose rows 
and columns are indexed by the elements G and whose 
entries are 

, a \ _ 1 if (g,h) e Ri 
if (g,h)^R i - 

So i 'th adjency matrix Ai is a 0, 1 matrix. 

Proposition 27 H3/ (i) A\ = 1, the identity ma- 
trix. 

fii) Sfe=l = J> where J is the matrix whose entries 
are all 1. 

( Hi ) A k = Ay for some k' G {1 , r} . 

(iv) AiAj = J2l=iP^k ViJ,k£ {l,...,r}. 

(v) p% = p*. Vi, j, fc e {1, r} = 
i^-Aj Vi,je{l,...,r}. 

The matrix M 15 may be written as a linear combina- 
tion of the adjency matrices in the following way 

Proposition 28 

r 

M D =J2x D (C i )A i . 

i=l 
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